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Abstract 

In the present work we examine both the hnear and nonhnear properties of two related 
PT-symmetric systems of the discrete nonhnear Schrodinger (dNLS) type. 

First, we examine the parameter range for which the finite PT-dNLS chains have real eigen- 
values and PT-symmetric linear eigcnstates. We develop a systematic way of analyzing the non- 
linear stationary states with the implicit function theorem at an analogue of the anti-continuum 
limit for the dNLS equation. 

Secondly, we consider the case when a finite PT-dNLS chain is embedded as a defect in the 
infinite dNLS lattice. We show that the stability intervals of the infinite PT-dNLS lattice are 
wider than in the case of a finite PT-dNLS chain. We also prove existence of localized stationary 
states (discrete solitons) in the analogue of the anti-continuum limit for the dNLS equation. 

Numerical computations illustrate the existence of nonlinear stationary states, as well as the 
stability and saddle-center bifurcations of discrete solitons. 

1 Introduction 

The subject of PT-symmetry and its physical implications has gained a tremendous momentum 
over the past few years. This field 'was initiated by the original proposal of C. Bender [10] who 
suggested that the linear Schrodinger operator with a complex-valued potential, which is symmetric 
with respect to combined parity (P) and time-reversal (T) transformations, is guaranteed to have 
real spectrum at a certain parametric regime. Thus, this was proposed as a viable alternative 
for the standard Hermitian quantum mechanics. Yet, it was the pioneering work in the group 
of D. Christodoulides both at the theoretical \12>\ [28] and experimental [29] levels that showcased 
nonlinear optics as a fertile ground for the physical implementation of the PT-symmetric potentials. 
These efforts have motivated a wealth of recent works, especially on the physical side, addressing 
various aspects of continuous and discrete PT-symmetric systems. These include among others the 
study of the fragility of PT-symmetry in linear problems [9l [26] , nonlinear stationary states of few 
site configurations (also referred to as oligomers, or plaquettes in two-dimensional lattices) [20 } [21 } 
[30l [32l [M] , as well as solitary waves and breathers in infinite systems both continuous [H [21 [51 [131 
[H [24] and discrete [121 [H [31] . 

While the number of studies of such PT-symmetric systems both in optics [18] and in atomic 
physics [151 [16] is rapidly growing, the volume of related mathematical works is rather limited and 



1 



mostly constrained to linear problems [H \7\ [22l [33] . It is the purpose of this paper to provide a 
number of rigorous results on nonlinear stationary states in PT-symmetric discrete systems. Our 
emphasis will be two-fold. 

First, we will consider finite PT-symmetric chains of the discrete nonlinear Schrodinger (dNLS) 
type [T7j . We examine their phase transitions (from a PT-symmetric oscillatory phase to the expo- 
nentially growing phase) when the gain and loss parameter is increased. The nonlinear stationary 
states bifurcate from the linear PT-symmetric states by means of a standard local bifurcation. 
On the other hand, we will also consider large-amplitude stationary states in an analogue of the 
well-known anti-continuum limit for the dNLS equation [25\ [27] , through a suitable rescaling of the 
PT-dNLS equation. This rescaling enables us to use the implicit function theorem to continue sta- 
tionary states from the limit, where they are effectively uncoupled and the gain and loss parameter 
is negligibly small. 

Second, we consider the case where the finite PT-symmetric chains are embedded in the infinite 
nonlinear lattice of the dNLS type. Again, we will examine phase transitions of such systems and 
will prove that the infinite PT-dNLS lattice has a wider stability interval compared to the isolated 
PT-dNLS chains. We also develop a proof of the existence of localized stationary states (discrete 
solitons) in the PT-dNLS equation. Numerical computations illustrate the theoretical results on 
existence of nonlinear stationary states, as well as the stability and saddle- node bifurcations of 
discrete solitons. 

Note that our technique allows us to prove existence of discrete solitons in the infinite PT-dNLS 
equation, but such discrete solitons are unstable because the phase transition in this infinite lattice 
occurs already at the zero value of the gain and loss parameter |26j . Earlier, existence of such 
discrete solitons was observed in numerical continuations from the diatomic PT-dNLS lattice |19] . 

The article is structured as follows. Section 2 covers fundamentals of the PT-dNLS equaiton. 
Section 3 is devoted to finite dNLS chains with four sections on eigenvalues of the linear PT-dNLS 
equation, local bifurcations of stationary states, bifurcations of large-amplitude stationary states, 
and numerical results. Section 4 is concerned with the PT-symmetric defects in infinite dNLS 
lattices and contains three sections on eigenvalues of the linear PT-dNLS equation, bifurcations 
of discrete solitons from the anti-continuum limit, and numerical results. Section 5 concludes the 
article with a summary and a discussion of future directions 

Acknowledgments: The authors thank James Dowdall for help at an early stage of the 
project during his NSERC USRA work and Dimitri Frantzeskakis for discussions on the subject 
of PT-symmetry. The work of P.K. is partially supported by the US National Science Foundation 
under grants NSF-DMS-0806762, and NSF-CMML1000337, from the Alexander von Humboldt 
Foundation and from the US AFOSR under grant FA9550-12- 1-0332. The work of D.P. is supported 
in part by NSERC and by the ministry of education and science of Russian Federation (Project 
14.B37.21.0868). 

2 Formalism of the PT-dNLS equation 

We consider the discrete nonlinear Schrodinger (dNLS) equation with non-conservative terms that 
introduce gains and losses of nonlinear oscillators. When gains and losses are combined in a 
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compensated network, the model referred to as the PT-dNLS equation takes the form 

= Un+l - '^Un + Un-l + i7(-l)"lin + |n„pU„, (1) 

where parameter 7 stands for the gain and loss coefficient. The finite PT-dNLS chain is defined for 
n S {1, 2, 2A^} for a positive integer N subject to Dirichlet boundary conditions uq = U2N+1 = 0, 
whereas the infinite PT-dNLS lattice is defined for all integers n on Z subject to the decay of ii„ 
to zero as \n\ — t- 00. The amplitudes n„ for all admissible values of n are complex- valued functions 
of time t. 

For notational consistency, we denote the sequence {un}nez of complex-valued amplitudes Un 
by the vector notation u. These vectors are considered in Hilbert space /^(Z) equipped with the 
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inner product (u, v) := "^nGZ^nVn and the induced norm ||u|| := (Xlnez 1""™!^) 
Let us formulate the evolution problem ([T|) in the complex Hamiltonian form 

i^ = Vflif(u,u), (2) 

where u is a collection of amplitudes Un for all admissible values of n (denoted by S) in some 
function space (denoted by X), the bar denotes complex conjugation, and the complex- valued 
Hamiltonian functional H : X ^ C takes the form 

H{U, U) = - ^ \Un+l - Un\^ + i^Y.{-lT\Un\^ + ^ '"""l^- 
neS n&S neS 

The dynamical system ([2|) is said to be PT-symmetric if there is a linear real- valued t-independent 
operator P : X ^ X such that 

p2 = / and N(u,u) = PN(Pu,Pu), (4) 

where N(u, u) := Vui^(u, u) and / : X — )■ X is an identity operator. 

If u(t) is a solution of the PT-symmetric dynamical system ^ for t in a symmetric interval 
J := {—to, to) C M with some positive to, then v(t) := Pu{—t) is another solution of the same 
system for t £ J. This statement can be checked by direct substitution. This symmetry suggests 
the following definition of the operator T : C{J, X) — ?• C( J, X): 

Tu{t) := u{-t), t G J. (5) 

Note that the operator T is sesquilinear in u and nonlocal in t. The letters P and T stand for 
parity and time reversal transformations, which correspond to fundamental symmetries in physics. 

When the vector field N is linear and given by N(u, u) := Tiu associated with a linear complex- 
valued bounded operator % : X ^ X, then the PT-symmetry is expressed in the standard form 

n = PUP. (6) 

Our first result is to show that the dNLS equation with compensated gain and loss terms ([T]) is 
a PT-symmetric dynamical system both for finite and infinite chains. 
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Lemma 1. Define S := {1, 2, 2A''} for a positive integer N and X := C^^. Then, the dynamical 
system ^ with complex-valued Hamiltonian |^ is PT-symmetric with respect to the operator P : 

[Pu]„ = U2iv+i-„, n G 5. 
Proof. We verify the statement with the exphcit computation. Given the definition of P, we obtain 

[P^u]„ = [Pu]2Ar+i_„ = Un 
2U2Ar+l-n + U2N-n + ^7(~ l)"^^2Ar+l-n + \u2N+l-n\'^U2N+l-n- 
= Un+1 - 2Un + Un-1 + Z7(- 1)^^"^^"'^U„ + \Un\'^Un 

= [N(u,u)]„, 

which recovers the second identity ([H). □ 

Remark 1. The symmetry of LemmaUl can be proven by simple reflection arguments. If the chain 
of oscillators has the damped site at the left end and the gained site at the right end, then since P 
reflects all oscillators about the middle point, the reflected chain has now the gained site at the left 
end and the damped site at the right end, that is, the reflected chain is equivalent to the complex 
conjugate chain. 

Remark 2. Although P in Lemma{l\ represents the fundamental physical symmetry, other choices 
of operator P are possible for the linear terms of the dynamical system ^-(3^. For instance, if 
N = 2, there exists another operator Pa such that P^ = id and H = PaTiPa, where 

-2a a 
-2a -a 
-a -2a 
a -2a 

with either a = -^ora = —-^ (this statement can be easily checked by means of symbolic software). 
Nevertheless, the operator Pa does not represent the PT-symmetry of the full nonlinear system 
|3i because the nonlinear term Nnon(u, u) := N(u, u) - Tiu does not satisfy the second identity 
For instance, we have 

[PNnon(Pu,Pu)]i = a'^ [2\2ui + U3\H2U1 + Us) + \ui - 2u3Wui - 2U3)] 

lr..o ii9 9— 9— ii9 ii9n 

= — |_17|iii| ui + 12|iiipu3 + 8W3W1 + GuiU's + 16|ti3| ui — 6\u^fUs\ 

/ Iml^Ul = [Nnon(u, u)]i, 

hence the second identity ^ is not satisfied. 



and 

[N(Pu,Pu)]„ = U2N+2-n " 

Applying P again, we obtain 

[PN(PU,PU)]„ : 



-2 - i-f 
1 





1 

-2 + ^7 
1 






1 

-2-i7 
1 





1 

-2 + ^7 



Pr, 
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Corollary 1. Let S = Z and X = l'^{Z,C)- For any fixed no G Z, the dynamical system ^ with 
complex-valued Hamiltonian ^ is PT-symmetric with respect to the operator P : Z^(Z) — t- /^(Z) 
given by 

[Pu\n = Uno-n, n £ Z. 

Proof. The proof follows from the proof of Lemma [1] when 5 = {1, 2, 2N} is replaced hy S = Z 
and the value of no G Z is arbitrary. □ 

3 Finite PT-dNLS lattices 

We shall now consider the PT-dNLS equation ([1]) for the finite chain Sn ■= {1,2, ...,2A^}, where 
is a positive integer, subject to the Dirichlet boundary conditions uq = U2n+i = 0. We study 
2N eigenvalues of the linear PT-dNLS equation to find the phase transition threshold 7Ar, which 
separates the neutral stability of the zero solution for 7 G {—'JNtIn) and the linear instability of 
the zero solution for I7I > 77V. We show that 7Ar is a monotonically decreasing sequence of such 
that 71 = 1 and — )■ as — )■ 00. 

We consider local bifurcations of nonlinear stationary states of the PT-dNLS equation ([T]) 
from the linear limit and prove that every simple eigenvalue of the linearized PT-dNLS equation 
generates a unique (up to a gauge transformation) family of the PT-symmetric stationary states 
in the parameter space. For 7 inside the stability interval (— 77v,7Ar), this yields the existence of 
2A^ branches of stationary states. These 2A^ branches are extended towards a large-amplitude limit 
with some intermediate bifurcations. 

We characterize the number and properties of the branches of the stationary states in the large- 
amplitude limit and show that there exist 2^ distinct branches for any 7 G (—71,71) = (—1, 1), for 
which |unp is large for all n G S^. We also discuss existence of other branches of the stationary 
states, which are centered at the middle sites of Sjy and for which \ui\'^ is small in the large- 
amplitude limit. 

These analytical results are illustrated with numerical approximations of the nonlinear station- 
ary states of the PT-dNLS equation ([1]) for the finite chain with N = 1,2, 3. 

3.1 Eigenvalues of the linear PT-dNLS equation 

We consider the linear stationary PT-dNLS equation on a finite chain Sn ■= {1, 2, 2A^}: 

EWn = Wn+l+Wn-l + i^{-lTWn, UGSn, (7) 

subject to the Dirichlet boundary conditions wq = W2N+1 = 0. Compared to the PT-dNLS equation 
([1]), the diagonal term of the discrete Laplacian operator has been included in the definition of the 
parameter E (see Remark [3] below). We shall find all 2A^ eigenvalues of the linear stationary dNLS 
equation ([7]) in explicit form, a result from which the phase transition threshold 7Ar is computed 
also explicitly. 
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Theorem 1. Eigenvalues of the linear eigenvalue problem ^ are found explicitly from the set of 
quadratic equations: 

^'^ + E^ = icos''(-^), l<j<N. (8) 



In particular, all eigenvalues are simple and real for 7 G (— 7Ar,7Ar), where 



7iv:=2cos(^^). (9) 



Proof By writing 

Xk = W2k-i, Vk = W2k, l< k< N, 
we can rewrite the linear eigenvalue problem ([Tj) in the equivalent form: 

Ex, = y,.,+y,-i^Xk, (10) 
Eyk = xk + Xk+i + tjyk, 

where the boundary conditions are now yQ = and x^r+i = 0. Expressing from the second 
equation of the system (jlOp and substituting it to the first equation of the system, we obtain a 
second-order difference equation 

(7^ + E'^)xk = Xk-i + 2xk + Xk+i, I <k < N, 

where the boundary conditions are now xq = —xi and xat+i = 0. Using the discrete Fourier 
transform, we represent the eigenvector satisfying the boundary condition xat+i = in the form 

Xk = sme{N+l-k), l<k<N. 

Parameter 6 in the fundamental interval [0, vr] defines uniquely the spectral parameter z := 7^ + 
from the dispersion relation 

Q 

z ■= j'^ + E"^ = 2 + 2cose = Acos^ -. (11) 

From the remaining boundary condition xq + xi = 0, we obtain 

6(1 + 2N) 
sm cos 2=0' 

where cos | 7^ (since {a^fclfcLi must not be identically zero). From the roots of sin ^(^+^^) ^ -^^g 
obtain the admissible values of 6 as follows: 

27r7 

l<J<iV, 



l + 2iV' 

which yields the result by ([TT]) . □ 
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Remark 3. For each eigenvalue E of the linear stationary dNLS equation ^ with the eigenvector 
w, there exists another eigenvalue E with the eigenvector P\v. This is an elementary conse- 
quence of the PT- symmetry, which produces a new solution v(t) = Pu(-t) = (Pw)e-*(^-2)* of the 
time- dependent dNLS equation (CP from the solution u(t) = we~*^^~^^* of the same equation. In 
particular, if E is a simple real eigenvalue (as in TheoremU^, then the eigenvector w can be chosen 
to satisfy the PT-symmetry 

W = Pw =^ Wn = W2N+l-n, n S Sn- (12) 



We list some numerical values of the phase transition thresholds: 



71 


= 2 cos 


vr 

3 ~ 


1, 


72 


= 2 cos 


2tt 
~5 " 


a 0.618, 


73 


= 2 cos 


Svr ^ 

T ' 


a 0.445. 



Note that lim7v-j.oo 7n = 0. 

3.2 Stationary states: local bifurcations 

We shall now consider nonlinear stationary states on a finite chain Sn, which satisfy the nonlinear 
stationary PT-dNLS equation: 

EWn = Wn+l+Wn-l+iji-l)"'Wn + \Wn\'^Wn, UGSn, (13) 

subject to the Dirichlet boundary conditions wq = W2N+1 = 0. We shall work in the space X = C^^. 

Assuming that the linear stationary PT-dNLS equation ([7]) admits a simple real eigenvalue 
Eq with the eigenvector wq £ X, we shall prove the existence of a branch of the PT-symmetric 
stationary states w G X satisfying the nonlinear stationary PT-dNLS equation (113p for £^ in a one- 
sided neighborhood of Eq. The solution branch is unique up to a gauge transformation: w — >■ e*"w, 
where a € M. This result corresponds to the standard local bifurcation of the nonlinear state w 
from the linear eigenstate wq, which is complicated here due to the presence of the PT-symmetry. 

The local bifurcation results were considered with formal perturbation expansions by Zezyulin 
&: Konotop p4j. Here we give a rigorous version of the same result. 

Theorem 2. Assume that Eq is a simple real eigenvalue of the linear stationary PT-dNLS equation 
Ij^ with the PT-symmetric eigenvector wq = Pwq in X = C^^. Then, there exists a unique (up 
to a gauge transformation) PT-symmetric solution w = Pw of the nonlinear stationary PT-dNLS 
equation US]) for real E > Eq. Moreover, the solution branch is parametrized by a small parameter 
a such that the map R 9 a — )• {E, w) gMx X is C°° and for sufficiently small a, there is a positive 
constant C such that 

W-Nf + \E - EQ\<Ca^. (14) 
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Proof. We write the nonlinear stationary PT-dNLS equation (jl3p in the abstract form 

(S-?^)w = Nnon(w), (15) 

where Ti : X ^ X is the hnear (matrix) operator associated with the right-hand side of the 
hnearized stationary PT-dNLS equation ([7]) and Nnon(w) : X — )• X is the cubic nonlinear part. 
We note that according to our assumptions, we have 

Ker(£'o — 'H) = span(wo), Ker(£^o — T~i-)^ = span(Pwo), 

where 

{Eq -n)+ = Eo-n = p{Eo - n)p. 

Using the standard Lyapunov-Schmidt method, we write 

E = Eq + ^, w = awo + u, (Pwo, u) = 0, (16) 

where (A, a, u) G C x C x X are determined from the nonlinear equations (llSp projected to Ker(£'o — 
H)^ and Ran(£'o — H)"*". Recall that by the Fredholm theory, Ker(£^o — "H)"*" is orthogonal to 
Ran(£;o - n) so that u G Ran(£;o - T-L). 

The projection to Ker(So ~ T~L)^ is written in the scalar form: 

Aa(Pwo,wo) = (Pwo,Nnon(awo + u)). (17) 

By the implicit function theorem, the projection to Ran(So ~ 'H)^ (not written here) guarantees 
the existence and uniqueness of a smooth (C°°) map from (A, a) £ to u G Ran(£'o ~ ^) C X. 
Moreover, for small values of A and a, there is a positive constant C such that 

||u|| <C(l + |A|)|a|^ (18) 

For a = 0, we have a unique zero solution u = and the equation (|17p is satisfied identically. 
In what follows, we assume a 7^ 0. 

We claim that (Pwq, wq) 7^ under the assumption that Eq is a simple eigenvalue of ^. Indeed, 
if (Pwo,wo) = 0, there exists a generalized eigenvector wi G X for the same eigenvalue Eq from a 
solution of the inhomogeneous equation 

(£'0 - T-L)mvi = -Wo, 

which is a contradiction to the assumption that Eq is a simple eigenvalue of T-L. 

Therefore, (Pwo,wo) ^ 0. Then, there exists a unique smooth map from a G C to A G C 
solving the bifurcation equation (117p . Moreover, for small values of a, there is a positive constant 
C such that 

|A(Pwo,wo) - |a|2(Pwo,N„on(wo))| < C\a\\ (19) 

Note that both (Pwo,wo) and (Pwq, Nnon(wo)) are real because of the PT symmetry of the 
eigenvector wq = Pwq and the nonlinear field Nnon satisfies the second identity Indeed, we 
have 

2N N 

{Pwo,^o) = (wo. Wo) = ^(wo)^ = [(wo)^ -I- (wo)^] 

n=l n=l 
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and 



2N N 
(Pwo,Nnon(wo)) = (wq, Nnon(wo)) = ^ |(wo)„P(wo)^ = ^ |(wo)nP [(wq)^ + (wq)^] . 

n=l 71=1 

Therefore, A is real at the leading order 0(1 ap). To exclude the gauge transformation, let us 
consider the real values of a. Because the nonlinear vector field Nnon preserves the PT-symmetry, 
the unique solution for u and A is PT-symmetric, so that u = Pu and A is real. The bound (fT4 
follows from (fT6]l . (fTSl) . and (fT9l) . To be precise, we obtain 

2 , /n/ 4^ A ._ (^'wo,Nnon(wo)) 



A = A2a' + 0(a^), As := 



(Pv^^o, Wo) 



It remains to prove that A2 > 0. However, using the explicit representation from Theorem [H for 
the eigenvalue with 6 = n ' 1 — i — obtain the eigenvector wq with components 

r- 27rj(iV + l-A;) ^- 27rj(N + 1/2 - k) 



Therefore, 



. 4 27ri(jV+l-fc) 



2 +7 ^ 2nj(N+l^k) 

Z^k=l ^^^^ l+2Af 



and the proof of the theorem is complete. □ 

Remark 4. The local bifurcation results do not apply in the limit N ^ 00 because of two reasons. 
First, the spectrum of the linear stationary dNLS equation ^ becomes continuous as N ^ 00. 
Second, for any 7 7^ 0, the spectrum includes complex (purely imaginary) points of E because 
7Ar — as N ^ CO. 

Let us consider the simplest example = 1 when Theorem [2] works. The two simple eigenvalues 
are Eq = ity^l — 7^ and the eigenvectors wq = -Pwq are given by the same expression 



il 

s/Eo + ^7 

In this case, A2 = |, so that we have the expansion 



wo = — 



E = ±^/l^ + A^ + 0{A^), A:=^a. 
In fact, it follows from the exact solution (|3U|) below that the error term 0{A'^) is identically zero. 
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3.3 Stationary states: bifurcation from infinity 

We shall now consider the stationary states of the nonlinear stationary PT-dNLS equation (jl3p in 
the limit of large values of E. This corresponds to the anti-continuum limit of weak couplings in 
the PT-dNLS lattice after a suitable scaling transformation (which is also discussed in [H]). Note 
that the standard anti-continuum limit arising when the coupling parameter in front of the discrete 
Laplacian operator vanishes fails to generate any solutions of the stationary dNLS equation (fT3]) 
for real values of E and 7 7^ 0. 

We shall develop methods to analyze a bifurcation from infinity for solution branches. In 
particular, we shall prove the existence of 2-^ branches of the PT-symmetric stationary states w of 
the nonlinear stationary PT-dNLS equation ()13p for 7 G (—1, 1) and for large values of E, for which 
\wn\'^ is large for all n £ Siy. The solution branches are unique up to the gauge transformation 
w — )• e*"w with a G M. The complication of proving this result is caused by the degeneracy of 
asymptotic solutions of the nonlinear algebraic system (jl3p in the limit E — )• 00. Indeed, setting 
w = \/^W and taking the limit — )• 00, we obtain an uncoupled set of algebraic equations with 
N PT-symmetric solutions 

where (/J^ E M is an arbitrary parameter and is a unit vector on the finite chain Sn. However, 
the space of solutions of the nonlinear algebraic system (jl3p in the limit E ^ 00 does not enjoy 
the linear superposition principle and parameters {^k}k=i ™ust be fixed from 0{l) conditions as 
£^ — 7- 00. To prove persistence of continuations of the limiting roots for large but finite values of E, 
we have to unfold the degeneracy of the nonlinear system by a special transformation, after which 
the result is guaranteed by the implicit function theorem. Along these lines, we prove the following 
main result. 

Theorem 3. For any 7 € (—1,1), the nonlinear stationary PT-dNLS equation il3\) in the limit of 
large real E admits 2^ PT-symmetric solutions w = Pw (unique up to a gauge transformation) 
such that, for sufficiently large E, the map E ^ w is C°° at each solution and there is a positive 
E-independent constant C such that 



^ - 2NE 



< C. (20) 



Proof. We set E = ^ and w = ^ for small positive 6 and write the stationary dNLS equation (jl3p 
in the equivalent form: 

{l-\Wn\^)Wn = 5{Wn+l + Wn-l+i7{-irWn), U^Sn, (21) 

subject to the Dirichlet boundary conditions Wq = W2N+1 = 0. We consider a PT-symmetric so- 
lution with W = PW such that the system can be closed at algebraic equations for 1 < n < 
subject to the reflection boundary condition W^^i = Wn. Note that parameter 7 G M is fixed. 

Case A^ = 1: In this case, we only have one nonlinear algebraic equation to solve: 

(1 - |VFip)VFi = 6[Wi- ijWi] . (22) 
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Setting Wi = 6**^^, we separate the real and imaginary parts of equation (j22p as follows: 

A\ = \ — (5cos(293i), — sin(299i) — 7 = 0. 

For any 7 G (—1,1), there exist two solutions for in [0, vr] from the second equation writ- 
ten as sin(2(^i) = —7. For each yp, we have a unique solution of the first equation written as 
= 1 =F b\J\ — 7^, from which we see that j4i = 1 + 0{f) as 5 — )■ 0. 

Case > 2: Let us now unfold the degeneracy of the nonlinear algebraic system (j2ip in the 
limit 5 — )• by using the transformation: 

, W3 = (^l^2^3)^/^e*^'+*^'+''^^ (23) 

, Wn = {A1A2 ■ ■ ■ ^^)i/2ei'/'i+i'P2+-+w^ 

where amplitudes Ai,A2,...,A]y and phases ifi, ip2, ^Pn are all real. After substitution and 
separation of real and imaginary parts, we obtain equations for phases 

7 = 0, 

— 1/2 

A^ sin(v92) + 7 = 0, 

A3 '/%in(<^3) - 7 = 0, (24) 

, - sin 2(991 + (^2 H h 1/9 at) - A'^^^'^ sm{ipN) + (-1)^7 = 0, 

and N equations for amplitudes 

/ 1 /2 

1 — Ai = 6A2 COs((/92), 

1 - A1A2 = 6{AI^^ C0S((/P3) + ^2 C0S(992)), 

< l-AiA2A3 = 6{Al^^cos{yD4) + A~^^\os{ip3)), (25) 

— 1/2 

^ 1 - ^1^2 • • • Aat = (5(C0S2((^1 + (/?2 H \- ^n) + Aj^ ' COs{ipN))- 

For (5 = 0, the system of amplitude equations ([25]) has a unique solution at the point Ai = A2 = 
■ ■ ■ = A]\[ = 1. The vector field of the nonlinear system is smooth with respect to {Ai,A2, ■ ■ ■ , A]\f) 
and S near this point for all {(pi,(p2, • • • , ^n) S T^, where T denotes the fundamental interval [0, 27r] 
subject to the periodic boundary conditions. The Jacobian matrix with respect to {Ai,A2, ■ ■ ■ , Ajsf) 
at this point has eigenvalue 1 of geometric multiplicity one and algebraic multiplicity N. By the 
Implicit Function Theorem, for all {(pi, ip2, ■ ■ ■ , ^Pn) £ and small 5 G M, there is a unique solution 
of the nonlinear system ([2^ such that the map {tpi, ip2, ■ ■ ■ , ^n,S) (^15^2! ■ ■ ■ , An) is and 
there is a positive 5-independent constant C such that 

|Ai - 1| + 1^2 - 1| + • • • + I^TV - 1| < C|<5|. (26) 



^2 sin(v92) - 
^3^^sin(v93) - 
< sin(</54) - 
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Bound (j20p follows from this bound and the scaling transformation. 

Now we consider the system of phase equations (I24p . which is 6 independent. Nevertheless, it 
depends on 5 via amplitudes {Ai,A2, . . . ,A]\[). For 5 = 0, the nonlinear system (|24p can be written 
in the explicit form: 



sin(992) = 7, 

sin(</53) = -7 + sin(v32) = 0, 
sin(v94) =7 + sin((^3) =j, 



sin 2(991 + (/J2 H 1- ^n) = (-1) 7 - sin{(p 



(27) 



Denote := 2((/Ji + (/J2 + • • • + ^n)- For any 7 G (—1, 1), there are 2^ possible solutions of (f27l) 
for {ijj,ip2, . . . ,^n) S T^, depending on the binary choice of the roots of the sinusoidal functions 
on the fundamental period. Because (fi = ^ — — • • • — , there are actually four solutions for 
ipi in T, however, the solutions with (pi S (vr, 27r] are reducible to the solutions with (pi E (0, vr] 
by the transformation W — t- — W, which is a particular case of the gauge transformation. In what 
follows, we only consider the two possible solutions for ipi in [0, vr]. 

The vector field of the nonlinear system ()24p with {Ai,A2, . . . ,A]^) obtained from the nonlinear 
system ([^^ is smooth with respect to ((/?i, ip2, ■ ■ ■ , ^pn) and 6. The Jacobian matrix with respect 
to (ifi, ip2, ■ ■ ■ ,(Pn) for 5 = is given by the matrix 









COS (932) 

-COs((/92) 







C0S((^3) 
-C0S(V?3) 









-2cos{ip) — 2cos(V') — 2cos(V') — 2cos(V') 










■C0S(V37V) 



Now it is clear that cos{ipn) 7^ for all 2 < n < if 7 G (—1, !)• In addition, the last equation in 
the system ()27p is given by either sm{ijj) = if is even or sm{ip) = 7 if is odd. In either case, 
cos{'tp) 7^ if 7 G (—1, 1). Hence, the Jacobian matrix is invertible if 7 G (—1, 1). By the Implicit 
Function Theorem, for all small 5 G M, there is a unique continuation of any of the 2^ possible 
solutions {Lp\,ip2, ■ ■ ■ jV^at) of the nonlinear system (|27p as a solution of the nonlinear system (I24p 
such that the map 6 — t- {(pi,Lp2, ■ ■ ■ ,(Pn) is C°° and there is a positive (5-independent constant C 
such that 

+ \^2-^*2\ + --- + \vn-^*n\<C\S\. (28) 



This completes the proof of the theorem. 



□ 



Remark 5. The number of solution branches grows as N ^ 00 for any fixed value of 7 in the 
interval (—1,1). However, all these solution branches are delocalized in the sense that |w„p w E 
as E ^ 00 for all n in Sn- Therefore, none of the solution branches of Theorem\M approach to a 
localized state (discrete soliton) as N ^ 00. 

Remark 6. Besides solution branches of Theorem\^ for any N >2 and 1 < M < N, there exist 
additional solution branches such that \wn\'^ k, E as E ^ 00 for N — M + \ < n < N + M and 
\wn\^ ~ as E ^ CO for 1 < n < N — M and N + M + l<n< 2N. These stationary states are 
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supported at 2M sites near the central sites in and their persistence is proved with a similar 
variant of the implicit function theorem (see the proof of Theorem\^ below). If N ^ oo, such 
stationary states approach to a localized state (discrete soliton). Note that the discrete solitons are 
unstable on the unbounded lattice because the continuous spectrum of the linearized dNLS equation 
^ is complex for any 7 7^ 0, recall that 77V — )• as N ^ 00. 

Remark 7. The arguments of the implicit function theorem can not be applied to construct solution 
branches which are centered anywhere but at the central sites in the set S^. Indeed, the numerical 
results below show that no such solution branches exist for large values of E. 



3.4 Numerical results 

We shall construct here the simplest nonlinear stationary states for N = 1,2,3. For = 1, this 
corresponds to the nonlinear dimer, where the solution branches can be obtained analytically, as 
in [201 [28l I32j . For N = 2, this corresponds to the nonlinear quadrimer and the solution branches 
can be at best approximated numerically [20^ I34j. For = 3, the numerical approximations of the 
nonlinear stationary states are added here for the first time. 

For = 1, we use the reduction W2 = wi and write wi = Ae~^'^ with real A and f. Then, the 
stationary PT-dNLS equation (fT3]) yields two equations 



sin(2(^) =7, = E- cos(2(^). (29) 
With two solutions of the first equation for 99 G [0, vr], we obtain two solution branches 



AI=Et^/T^, 7 g (-71,71), (30) 

where 71 = 1. The two solution branches coalesce into one branch for 7 = 71 and disappear via a 
saddle-center bifurcation for 7 > 71. 

Positivity of A"^ shows that E > E± = ity^ 1 — 7^, where E± are the simple eigenvalues of the 
linear stationary PT-dNLS equation for 7 G (—71,71)- We note that — t- as — t- E± and that 
A'^ ~ as — )• 00. These analytical results clearly illustrate the bifurcation results in Theorems 
El and El 

For N = 2,we use the reduction W4 = wi and = W2- Writing wi = Ae"**^"*^ and W2 = Be~^^ 
with real A, and in the nonlinear stationary dNLS equation (|13p . we obtain the system of 

nonlinear equations: 

sin((^) = ^, sin(2^) = i^-^2 (^1) 

and 

Ac" = AE - Bzo^{ip), = BE - Acos{ip) - Bcos{2iIj). (32) 



The roots of the algebraic system (I3ip and p2|) can be investigated numerically and the results 
depend on the value of 7. Figure 1 shows the solution branches on the {E, ^^)-plane (top) and the 
{E,B^)-plane (bottom) for 7 = 0.5 < 72 w 0.618 (left), 7 = 0.75 (middle), and 7 = 1.1. Note that 
no solution branches exist for 7 > 72 := 2 cos ^ ~ 1.618, because no simple real eigenvalues occur 
in the linearized dNLS equation ([7]) for these values of 7. 

According to Theorem [21 we count exactly four (2A^ = 4) solution branches for small amplitudes 
A and B for 7 < 72 and exactly two small solution branches for 72 < 7 < 72. According to 
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Figure 1: Nonlinear stationary states for N = 2 and for 7 = 0.5 (left), 7 = 0.75 (middle), and 
7 = 1.1 (right). The top and bottom rows show components \wi\^ and |w2p respectively. 




Figure 2: Nonlinear stationary states for = 3, and for 7 = 0.25 (top) and 7 = 1.1 (bottom). The 
left, middle, and right columns show components |ii'2pi and |w3p respectively. 
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Theorem [3l we count exactly four (2-^^^ = 4) solution branches for large amplitudes A and B if 
7 < 1, whereas all solution branches terminate before reaching large amplitudes if 7 > 1. 

The two branches for small values of A and large values of E are attributed to the solutions in 
Remark [6l The corresponding values of B are large. On the other hand, no branches exist for large 
A and small B as E gets large, see Remark [71 

For = 3, we write wi = Ae-'^'^+'l'+^^ = wq, W2 = Se-^W'+^) = W5, and W3 = Ce~'^ = m- 
The roots of the resulting algebraic system are investigated numerically by a homotopy method 
and the results are shown on Figure [2j We count six (2A^ = 6) branches in the small-amplitude 
limit if 7 < 73 ~ 0.445 and four branches if 73 < 7 < 73; = 2 cos ^ ~ 1.247. We also count eight 
^2Af=3 _ g-j branches for large amplitudes if 7 < 1 and no branches for large amplitudes A'^ 
if 7 > 1. More branches are counted for large amplitudes B^ and even more branches for large 
amplitudes C^. Overall, the results for iV = 3 are similar to the results for N = 2. 

4 PT-symmetric defects embedded in infinite PT-dNLS lattices 

We shall now consider an infinite PT-dNLS lattice, where the particular emphasis is on the existence 
and stability of localized stationary states (discrete solitons). Because the phase transition in the 
PT-dNLS equation ([T|) on the infinite lattice occurs already at jn^oo = 0, there is no way to obtain 
stable discrete solitons in such systems with extended gain and loss [26]. Therefore, we modify the 
PT-dNLS lattice by considering the PT-symmetric potential as a finite-size defect. Such defects 
were considered recently in the physical literature [31] and [H [5] . 

Let be a positive integer and Sn '■= {1, 2, 2iV} be the sites of the lattice, where the 
PT-symmetric defects are placed. The model takes the form 

^~lt' ^ ^""^^ ~ """"^ n(-l)"Xne5jv'"n + \Un\'^Un, n £ Z, (33) 

where XtigSn ^ characteristic function for the set Sn- When A'' = 1, the PT-dNLS equation (p3]) 
corresponds to the embedded dimer in the infinite PT-dNLS lattice. When = 2, it corresponds 
to the embedded quadrimer, and so on. 

We study the linearized PT-dNLS equation and find the phase transition threshold 77V. It is 
quite remarkable that 77V > Jn for any A^ G N, in particular, 71 = \/2. Nevertheless, is still a 
monotonically decreasing sequence of A^ such that 77V — )• as A^ — )■ 00. 

Then, we employ the large-amplitude (anti-continuum) limit of the PT-dNLS equation (j33p 
to study the existence of discrete solitons supported at the PT-symmetric defect Sm- For recent 
results on existence of discrete solitons in the anti-continuum limit for the regular dNLS equation 
(in the absence of PT-symmetry) , see e.g. [3l [11]. We find that for all 7 G (—1, 1), 2^ branches of 
the discrete solitons exist in this limit, for which is large for all n G Sn- 

The existence and stability of discrete solitons is illustrated numerically and we show that 
the stable branches of the discrete solitons for 7 7^ originate from the stable branches in the 
Hamiltonian version (7 = 0) of the dNLS equation \25\ 127] . 
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4.1 Eigenvalues of the linear PT-dNLS equation 

We consider the linear stationary PT-dNLS equation: 

EWn = Wn+l + Wn~l + i^i-l)"'XneSN'^n, n £ Z. (34) 

Because the PT-symmetric potential is compact, the continuous spectrum of the linear PT-dNLS 
equation (pl|) is located for E S [—2,2]. Besides the continuous spectrum, isolated eigenvalues 
may exist outside the continuous spectrum. To characterize isolated eigenvalues, we introduce a 
parametrization 

E:=2cos9, Re(6l) G [-tt, vr], Im(6') > 0, (35) 



and look for exponentially decaying solutions of the linear PT-dNLS equation (j34p in the form: 



"^n-i ^.._^_je{n-2N) n > 2N , ^^^^ 



which still leaves a set of 2N unknown variables {ifnlne^jv- To find {wn}n(^SNi close the linear 
eigenvalue problem at the algebraic system 

2 cos OWn = Wn+l + Wn-1 + i^{-lTWn, U £ Sn, (37) 

subject to the boundary conditions 

Wo = Wie'^, W2N+1 = W2Ne'^. 

Note that each eigenvalue E is complex if lm{6) > and there exists a complex conjugate 
eigenvalue E by the PT-symmetry (see also Remark [3]) . The following result is similar to the result 
of Theorem [H 

Theorem 4. A new symmetric pair of complex-conjugate eigenvalues of the linear PT-dNLS equa- 
tion [34^ bifurcates at \'y\ = ^N,k, where 

■ji(2k — 1) 

-fN,k ■■= 2 COS ^4^^, 1 < A: < iV, (38) 

and persists for \ j\ > 7jv,fc except possibly finitely many points on any compact interval of 'j, where 
the pair coalesces into a double (semi-simple) pair of real eigenvalues. In particular, no complex 
eigenvalues exist for ^ £ (—JNt^n), where 

7r(2A^-l) , , 

7JV := lN,k = 2 cos — . (39) 



Proof. We set 

Xk = W2k-i, Vk = W2k, I <k < N 
and rewrite the linear eigenvalue problem (|37p in the equivalent form: 

- 2cosex, = y,^, + y,-^JX,, ^^^^^ ^40^ 
2 cos dyk = Xk + Xk+i + tjyk , 
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subject to the modified boundary conditions uq = xie*^ and xat+i = yN^-^^ • Expressing from the 
second equation of the system (HOl) by 

Xk + Xk+l , ^ u ^ AT ^ 

yk = 7, a l<A;<iV-l, VN = —n -, 

2 cos a — I'y e — Z7 

and substituting these expressions to the first equation of the system, we obtain a second-order 
difference equation 

(7^ + 4 cos^ e)xk = Xk-i + 2xk + Xk+i, l<k<N, (41) 
where the boundary conditions are now 

xo = {e'' -ij)e''xu x;v+i = ^^. (42) 

e — ij 

Note that equations (j4ip are obtained by multiplying every term by 2 cos — i-f, which is hence 
supposed to be non-zero. 

The second-order difference equation (|4ip admits an exact solution 

Xk = C+e^"(*^-i) + C_e-^"('=-^) , (43) 

where a is defined from the transcendental equation 

72 + 4cos2 6l = 2 + 2cosa, (44) 

and (C+,C_) are non-zero solutions of the linear system following from the boundary conditions 

e-^"C7+ + e^"C_ = {e'^ - n)e'^ {C+ + C_), 
{e~'^ - i7)(C+e^"^ + C-c"^"^) = e^^(C+e*"(^-^) + C_e-*"(^"^'). 

Note that for fixed 7, the values of 9 are obtained from the characteristic equation for this linear 
homogeneous system, after the values of a are excluded from the transcendental equation (|44p . 
Also note that only the values of 9 with Im(^) > determine isolated eigenvalues of the linear 
stationary PT-dNLS equation by means of the representation ([35|) . 

After some algebraic manipulations, the characteristic equation for the linear system takes the 
form 

(cos(aA^) sin 9 sin a — i sin(aA^) cos 9{1 — cos a)) (7 + 2i cos 9) = 0. 

The equation 7 + 2i cos 9 = gives an artificial root corresponding to the value a = vr, because the 
second-order difference equation (|4ip is obtained by multiplying every term by 2 cos 9 — ij. There- 
fore, we drop this nonzero factor from the characteristic equation and reduce it to the transcendental 
equation 

^2ie ^ sin(iV + l)a-siniVa ^ cos (iV + ^) a 
^ sinA^Q - sin(iV - l)a cos {N - ^) a 
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Case A'' = 1: Equation (jl5]) yields e^*^ = 2 cos a — 1. When this constraint is used in equation 
(I44p . we obtain e~^*^ = 1 — 7^. Since 7 G M, we obtain the existence of a simple eigenvalue with 
Ke{6) > for 7^ > 2 and the bifurcation occurs for 7 = \/2 and corresponds to 6 = it^, when 

cos e = o. 



Case N > 2: In a general case, we first consider bifurcations of complex values of 9 from real 
values of 9. Hence, we set € R and realize from (j44p that Im(cosa) = 0, which implies either 
a G M or a G vrfc + zM for any k £ Z. In both cases, the characteristic equation ()45p with real 6 
implies that either 9 = 01 6 = ziz^. 

If 9 = 0, then the characteristic equation MSh reduces to the equation sin(aA^) = 0, whereas 
equation ()44p implies that 7^ = —2(1 — cos a) < 0, which is outside of the parameter range we 
are interested in. (Recall here that no eigenvalues with lm9 > exist in the self-adjoint case with 
7 = 0.) 

On the other hand, if 9 = it^, then the characteristic equation (j45p reduces to the equation 
cos(aA^) = (recall here that the artificial root a = vr is neglected) or equivalently, 

TT(2k-l) 

a = ak:= , 1< A; < iV. 

From equation ()44p . we obtain that the bifurcation occurs at 

'7r(2/c - ] 



7 



4 cos 



2 



4 cos 



2N 



which corresponds to the values (j38|) . 

Next, we show that the values with lm{9) > correspond to the range 7^ > 7|- This implies 



that a new isolated eigenvalue E of the linear stationary PT-dNLS equation (j34p with Im(£') ^ 
bifurcates from the value E = that corresponds to 9 = ±^ at 7 = ±7^ and persists for all I7I > 7^. 
A symmetric complex-conjugate eigenvalue E exists by the PT symmetry of the stationary PT- 
dNLS equation (IM]) . 



To show the above claim, we use the parametrization z := e*", C := e^*^, and s := 7^. The 



system of transcendental equations (j44p and (I45p becomes the system of algebraic equations: 



(2(1 + z 



2N-l\ ^2N+l 



Z 



s, 

1. 



We know that for s = sq := 7^, there exists a solution of this algebraic system for = —1 and 
-1. Therefore, we consider the continuation of this complex-valued root in real values of s. 



^2N 



By computing the derivative in s, we obtain 

-2 



2(1 + ^2^-1) 



1 - z- 

C(l + 2iV^2iV-i) _ 



{2N + l)z^^ 

For s = Sq, we obtain from this linear system that 

2Nz 





r K 1 




' 1 ' 




ds 






dz 









- ds - 





dC 

ds 



(i-zy 



N 



2sin^^) 
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and since C = e^*^) this proves that ■^lm(6)\g^j-n > 0. By continuity of roots of the algebraic 
system above, Im(^) remains positive for s > sq, that is, for 7^ > 7^, near the bifurcation point. 

We have shown that roots with lm{9) > bifurcate from the points 6 = ziz^ and remain in the 
upper half-plane for all 7^ > 7fc- These roots correspond to complex eigenvalues E if Re(0) 7^ (in 
which case E > 2) or if Re(0) 7^ ibvr (in which case E < 2). Both situations can not be a priori 
excluded, however, we have here two facts: 

• When Re(^) = or Re(0) = ivr, a pair of complex-conjugate eigenvalues E coalesce at the 
real line into a double (semi-simple) eigenvalue, because of the PT-symmetry with lm{6) ^ 
generates two linearly independent eigenvectors for the same real eigenvalue. 

• Roots 9 are analytic with respect to parameter 7 by analytic dependence of roots of the 
algebraic system. 

Combining these two facts together, we realize that that the double (semi-simple) eigenvalues E 
can not split along the real axis (as each real eigenvalue after splitting would then become a double 
eigenvalue by the PT-symmetry). And they can not persist on the real axis as a double eigenvalues 
because of analyticity of the parameter continuation of roots 9 with respect to 7. Therefore, these 
double real roots split back to the complex domain. In addition, analyticity of the parameter 
continuations guarantees that there are finitely many points on any compact interval in 7, where 
the pairs of complex-conjugate eigenvalues can coalesce at the real axis. The proof of the theorem 
is complete. □ 

We list some numerical values of the phase transition thresholds: 



71 


= 2 cos 


TT 

4 ~ 


^/2, 






Svr ^ 




72 


= 2 cos 




a 0.765 






Svr 




73 


= 2 cos 


12 ^ 


» 0.518, 



with lim7v-5.oo In = 0. Note that 77V > 77v for any N £N. 

4.2 Stationary states: bifurcations from the anti-continuum limit 

We now consider the stationary states, which satisfy the nonlinear stationary PT-dNLS equation 

EWn=Wn+l+Wn~l + h{-'^)"'Xn(^SNl^n + \Wn\'^Wn, n £ Z. (46) 

We explore the large-amplitude limit similarly to Section 3.3. Hence, we set E = and w = 
where (5 is a small positive number. The stationary PT-dNLS equation ()46p is rewritten in the 
equivalent form: 

{l-\Wn\^)Wn = 6{Wn+l + Wn-l + h{-irXneSMWn), n G Z. (47) 

We consider the PT-symmetric solutions with W = PW, where P is given by [PW]„ = W2N+i-n- 
Note that the choice of no in P given by Corollary [1] is adjusted to the center of the PT-symmetric 
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defect. Therefore, the existence of the PT-symmetric solutions can be considered in the framework 
of the following two subsystems: 

{l-\Wn\'')Wn = 6{Wn+l + Wn-l+il{-lTWn), l<n<N, (48) 

and 

(1 - \Wn\'')Wn = 5{Wn+l + t^n-l), « < 0, (49) 

subject to the boundary condition Wn+i = Wn- The following result is similar to the result of 
Theorem [3l 

Theorem 5. For any 7 G (—1,1), the nonlinear stationary PT-dNLS equation in the limit 
of small positive S admits 2^ PT-symmetric solutions W = PW G /^(Z) (unique up to a gauge 
transformation) such that, for sufficiently small 5, the map 6 — )• W is C°° at each solution and 
there is a positive 6-independent constant C such that 

|W - Wol < C6, (50) 

where Wq is a solution of Theorem\^ (after rescaling). 

Proof. We consider small solutions of the subsystem ()49p for a given Wi G C and small 5 G M. 
Parameter 7 G M is fixed. The nonlinear system represents a bounded map from (W_ , Wi, 6) G 
X C xR to where Z_ is the set of negative integers including zero. For 6 = and 

arbitrary Wi G C, W_ = is a root of the nonlinear map and the Jacobian with respect to W_ at 
W_ = is invertible. By the Implicit Function Theorem, for all Wi G C and small (5 G M, there is 
a unique solution of the nonlinear system (j49p such that the map {Wi,5) — t- W„ is and there 
is a positive (5-independent constant C such that ||W_||;2(^_) < C\S\. 

Substituting Wq from the map constructed above to the first equation of the subsystem (jl8|) . 
we close the system at N nonlinear equations for {Wn}i<n<N- The only difference from the 
nonlinear equations considered in the proof of Theorem [3] is the boundary condition for given Wq, 
however, Wq is small as Wq = 0{5). The two applications of the Implicit Function Theorem 
developed in the proof of Theorem [3] apply directly to our case and yield the assertion of this 
theorem. The bound (f50]) follows from bounds (i26]l and (1281) . □ 

Remark 8. A remark similar to Remark\^ applies on the infinite lattice as well. Besides localized 
states of Theorem\^ for any N >2 and \ < M < N , there exist additional soliton states such that 
\Wn\'^ 1 as 5 ^ for N - M + 1 <n < N + M and ^ as 5 ^ for n < N - M and 

n > N + M + 1. These stationary states are supported at 2M sites near the central sites in the set 
Sn- 

We give details of the perturbative expansions for the two (most fundamental) discrete solitons 
supported by the dimer defect for = 1. The subsystems (j48p and (j49p are rewritten explicitly as 
follows: 

{l-\Wi\'^)Wi = 6 {Wq + Wi - i-fWi) , (51) 

{l-\Wn\^)Wn = 6{Wn+l+Wn-l), n < 0. (52) 
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The perturbation expansion 



(53) 



allows us to compute at the leading order Wn = e where 9 G [0, vr] is arbitrary at this 

point. At the 0{6) order, we obtain the equations: 

W^'^ = W!^l + W!^l,, n<0 (54) 
- (e-2^^wf ^ + M^^')) = e^'-ne-^', n = l (55) 

Set Wi^'^ := Ze~^^ and rewrite (\55h as —{Z + Z) = e'^^^ — ij. The solvability condition is sin(20) = 7 
and it gives exactly two values for 9 € [0,7r] for any 7 G (—1, !)• Then, the first-order correction 
term is found explicitly as follows: 

^i'^ = -^e-'' cosi29)5n,i + e-''5n,2. 

The perturbation expansion (j53p can be continued to higher orders of 6 thanks to C°° smoothness 
in Theorem [51 Hence, we obtain two branches of soliton states supported by the dimer defect. 

Remark 9. The phase transition threshold for Si is 71 = ^/2, whereas the soliton states of Theorem 
are only constructed for 7 G (—1, 1) in the limit E — )■ 00. Numerical studies (see Figurel^ show 
that the solution states exist for 7 G ( — 1,1) for any fixed value of E. 

4.3 Numerical results 

We now test these analytical results for the PT-symmetric chain with embedded defects. 

First, we consider the linear limit of such chains and examine the corresponding PT-symmetry 
phase transitions expected to occur at ^m- In particular, in Fig. [3l we consider the case of = 1, 
i.e., a single embedded dimer which has been predicted to have a PT-phase transition at 71 = ^/2, for 
the infinite lattice. In the figure, we can clearly discern the relevant transition (the corresponding 
vertical dashed line shows the theoretical prediction of 71). Nevertheless, for the finite lattice 
considered (here 800 sites are used), an additional bifurcation is observed at 7 = 1 (see also the 
works of [5l [26]). This bifurcation is suppressed at the infinite lattice limit, but for a finite chain 
a "bubble" of complex eigenvalues arises around E = {) (which is the middle point of the spectral 
band). This bubble keeps expanding and slowly increasing in imaginary part between 7 = 1 and 
7 = ^/2 (notice that this growth is barely visible in the linear scale of the figure but it is noticeable 
in the logarithmic scale of the inset). At the latter critical point, the rapid growth of the isolated 
unstable eigenvalue pair becomes dominant for the instability of the lattice with an embedded 
dimer. 

Similar conclusions can be drawn for the case of embedded quadrimer (A = 2) and hexamer 
(A = 3) from Fig. [H The figures reveal, however, that in this case in addition to the actual (infinite 
chain) critical points of 72 ~ 0.765 and 73 ~ 0.518, respectively, there are multiple additional 
points of weak instability emergence due to finite size effects. Such features are noticeable due to 
bifurcations of instability bubbles at the edges of the spectral band (at = 2 and E = —2), at 
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Figure 3: The left panel shows the dependence of the maximal imaginary part Aj of the eigenvalues 
A = Ar + i\i versus 7 for iV = 1. While this is shown in a linear scale with the green dashed 
line denoting the bifurcation point 71 = y^, the inset in a semi-logarithmic scale demonstrates the 
finite-size instability emerging at 7 = 1, which is denoted by the red dashed line. The right panel 
shows the spectrum of the linear PT-dNLS equation (j34p for 7 = 1.43 (top) and 7 = 1.4 (bottom). 




Figure 4: The case of = 2 (left panel) and = 3 (right panel) for a finite lattice of 800 sites. 
The panels are similar to the left panel of Fig. [31 but feature two weak instabilities due to finite size 
effects (green and red dashed line), in addition to the strong instability at 7jv (blue dashed line). 
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7 ~ 0.46 for N = 2 and at 7 ~ 0.30 for A'^ = 3. Additional bubbles emerge near the middle point 
of the spectral band (at £; = 0), at 7 0.61 for iV = 2 and at 7 0.45 for = 3. 

Next, we turn to the existence of nonlinear stationary solutions in the PT-symmetric dimer 
(N = 1) embedded in the nonlinear chain, as shown in Fig. [5l We have identified two branches 
of localized states starting from the Hamiltonian limit where 29 = and 26 = tt (i.e., in-phase 
and anti-phase discrete solitons, respectively). Indeed, the left panel shows the phase difference 
between w\ and W2 for the PT-symmetric embedded dimer. These results clearly illustrate that 
at the Hamiltonian limit of 7 = 0, the system starts from the two well-known in phase and out- 
of-phase solutions [25]. The former (lower norm one) is unstable for our focusing nonlinearity, 
while the latter one is spectrally stable. Interestingly, in accordance with what is known for an 
isolated dimer [20^ [28l [3T] , these two solutions disappear in a saddle-center bifurcation at 7 = 1 . 
In fact, both this increased proximity and the eventual collision and disappearance are captured 
very accurately by the solvability condition sin(2^) = 7. The resulting angle from the numerical 
computation and from the essentially coincident analytical prediction are shown in the left panel 
of Fig. [5l Note that while the results are shown in Fig. [5] for 5 = 0.01, they remain similar not 
only for smaller values of 5 (such as e.g. 0.001), but even for larger values up to 5 = 1 that were 
considered. In particular, the point of the saddle-center bifurcation 7 = 1 has been found to be 
identical for other values of 5. This confirms the observation made in Remark [9l 

In the Hamiltonian case of 7 = 0, our linear stability results for these branches fall back on 
the analysis of [25]. In fact, we retrieve the exact same condition, for the leading order eigenvalue 
correction, namely /5 = 4cos(2^). This expression reveals the instability of the in-phase mode 
and the stability (for our focusing nonlinearity) of the out-of-phase one. In the presence of gain/loss 
(i.e., for finite 7), the fundamental difference lies in the existence condition which mandates that 
sin(20) = 7 and therefore the eigenvalues of the in-phase unstable state approach the origin (from 
the real axis), as do the ones of the out-of-phase marginally stable state (from the side of the 
imaginary axis) as 7 is increased towards 1. These eigenvalue pairs end up colliding at the origin 
at the point of the PT-phase transition at 7 = 1. The trajectories of the two sets of eigenvalues are 
shown in the right panel of Fig. [5l It can be seen that similarly to the Hamiltonian case of [25], 
the agreement is better for the unstable branch of real eigenvalues. Nevertheless, for both branches 
the comparison of computation and analysis is highly favorable. 

For the case of = 1, we now turn to numerical simulations, in order to briefly discuss the 
difference between the manifestation of the instability of the in-phase localized states for different 
values of 7. Two prototypical examples are shown in Fig. [6l Both panels illustrate the evolution 
of the two most central, maximum amplitude sites of the solution. In the Hamiltonian case of 
7 = shown in the left panel, it can be seen that small perturbations give rise to an amplified 
symmetry breaking in the dynamics. Nevertheless, the conservative nature of the ensuing dynamics 
ascertains the rapid saturation of this symmetry breaking and the eventual oscillations that arise 
lead to a (nearly) periodic alternation between a symmetric and a symmetry broken state. On the 
other hand, for 7 = 0.5 we observe a drastically different manifestation of this instability shown 
in the right panel. More specifically, the site associated with gain grows indefinitely in a nearly 
exponential form, as illustrated in the inset. At the same time, the site associated with damping 
decreases in amplitude in a similar fashion. 

Lastly, we address the bifurcation and stability results for the case of = 2 i.e., for an embedded 
quadrimer. In this case, the eigenvalues of the linearization problem are shown in Fig. [71 We have 
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Figure 5: The left panel shows the relative phase 26 between the two central sites obtained by the 
numerical computation and the solvability condition sm{29) = 7 (a green dash-dotted line). The 
right panel shows the squared eigenvalue of the linearized PT-dNLS equation at the discrete soliton 
versus 7 for each of the two branches. The stable (blue solid) branch of the anti-phase solution (at 
7 = 0) has a negative A^, while the unstable (red dashed) branch of the in-phase solution (at 7 = 0) 
has a positive A^. The corresponding theoretical predictions are shown by green dash-dotted lines 
in very good agreement with the numerical results. 
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Figure 6: The left panel shows the evolution of the unstable in-phase solution for 7 = 0. The 
amplitude of the central-most two nodes is shown as a function of time. The symmetry-breaking 
results from the amplification of a noise in the initial data and manifests the instability via a 
(nearly) periodic recurrence between symmetry-restored and symmetry-broken phases. The right 
panel shows the same instability for 7 = 0.5. The instability leads to an amplification of the gained 
site amplitude and a decay of the damped site amplitude. 
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Figure 7: Stability analysis results for the case of = 2. The left panel depicts the + + ++ solution 

(which has three real eigenvalue pairs) by solid (blue) lines, as well as the H h solution (which 

has two imaginary and one real pair) by dashed (red) lines. The right panel depicts the H 1 — 

solution (which has three imaginary eigenvalue pairs) by solid (blue) lines, as well as the + H 

solution (which has two real pairs and one imaginary) by dashed (red) lines. Each pair of these 
configurations collides and disappears at 7 = 1. 

examined four principal configurations with all four sites excited (there exist also configurations 
with two sites excited as for = 1, in accordance with Remark [8]). In agreement with Theorem [3l 

these four configurations are coded as + + ++ (when all sites excited in phase), H 1 — (when all 

sites are out of phase with their immediate neighbors), + H , and H h. As can be seen in 

the figure and as is known from the Hamiltonian limit of 7 = 0, the only spectrally stable among 

these configurations is the out of phase configuration H 1 — with three imaginary eigenvalue 

pairs, while the in-phase configuration + + ++ is the most unstable with three real pairs. Upon 
continuation over 7, for both of these configurations, two of the eigenvalue pairs move towards zero 
(which they reach as 7 — )■ 1), while one remains real for the in-phase, and imaginary for the out of 
phase. Remarkably, pairwise these configurations collide and disappear in the limit of 7 — ?• 1. More 

specifically, the in-phase configuration + + ++ collides with the configuration H h- Similarly 

the out-of-phase configuration H h— collides with the configuration + H . It should also be 

noted that we weren't able to continue any asymmetric mixed phase configurations (such as e.g. 

+ + H — , H , — I , 1 — or their opposite parity variants) past the Hamiltonian limit 

of 7 = 0. 

5 Conclusion 

In this paper, we examined two distinct scenaria for PT-symmetric dynamical lattices. In the first 
one, A^-site PT-symmetric chains were considered as a finite-dimensional dynamical system. In the 
second one, we have considered the embedding of the finite PT-symmetric system as a defect in an 
infinite dNLS lattice. In both cases, we have examined the linear problem, explicitly computing the 
corresponding eigenvalues and identifying the strengths of 7 beyond which instabilities (and the 
phase transition breaking the PT-symmetry) arise due to real eigenvalues. We have also considered 
the nonlinear states when they emerge from the linear limit, as well as when they arise from a highly 
nonlinear limit under suitable rescaling (analogous to the anti-continuum limit of the standard dNLS 
lattice). In that case, we argued about the disparity of the branch counts in these two limits (for 
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general N), which suggests the existence of a number of bifurcations, such as saddle-center ones, 
at intermediate values of the corresponding parameter. 

In the case of the infinite PT-symmetric PT-dNLS equation, 

=Un+l-2Un +Un-l+i'y {-!)"■ Un + \Un\'^Un, n G Z, (56) 

we note that the phase transition threshold is now set at 'Jn^oo = 0. Therefore, for any 7 7^ 0, the 
linear PT-dNLS equation is unstable with a complex-valued continuous spectrum. Nevertheless, 
we can still obtain existence of stationary localized states (discrete solitons) in the large-amplitude 
limit for 7 G (—1, 1) for any of the configurations described in Theorems [3] and [5l Moreover, the 
discrete soliton can be centered at any site uq £ Z because of the shift invariance of the PT-dNLS 
equation ([56|) . 

There are other numerous directions in which one can envision generalizations of the present 
study. In the present work, we considered the case where there is a single parameter 7, for each of 
the N pairs of sites with gain and loss. However, it is also relevant to generalize such considera- 
tions to the case of many independent parameters for such sites |34j . On the other hand, one can 
consider generalizations of the present setting that aim towards the case of higher dimensionality. 
Arguably, the simplest such generalization concerns the setting of two one-dimensional coupled 
(across each of their sites) chains in the form of a railroad track as in j30] and the consideration of 
multi-site excitations therein. However, the genuinely higher dimensional problem and the exam- 
ination of generalizations of plaquette configurations j21] . whereby the potential of vortices exists 
in the Hamiltonian limit is of particular interest in its own right. These themes will await further 
consideration. 

References 

[1] F.Kh. Abdullaev, Y.V. Kartashov, V.V. Konotop, and D.A. Zezyulin, "Solitons in PT- 
symmetric nonlinear lattices", Phys. Rev. A 83 (2011), 041805(R) (4 pages). 

[2] N. V. Alexeeva, I. V. Barashenkov, A. A. Sukhorukov, and Yu.S. Kivshar, "Optical solitons 
in PT-symmetric nonlinear couplers with gain and loss", Phys. Rev. A 85 (2012), 063837 (13 
pages) 

[3] G.L. Alfimov, V.A. Brazhnyi, and V.V. Konotop, "On classification of intrinsic localized modes 
for the discrete nonlinear Schrodinger equation", Physica D 194 (2004), 127-150. 

[4] J. D'Ambroise, P. G. Kevrekidis, and S. Lepri, "Asymmetric wave propagation through non- 
linear PT-symmetric oligomers", J. Phys. A Math. Theor. 45 (2012), 444012 (16 pages) 

[5] J. D'Ambroise, P.G. Kevrekidis, and S. Lepri, "Eigenstates of lattices with embedded PT- 
symmetric or Hamiltonian oligomer defects" , I a r Xiy : 1 2 1 1 . 5 70 7l . 

[6] T. Azizov and C. Trunk, "On domains of PT symmetric operators related to —y"{x) + 
(-l)"x2"2/(x)", J. Phys. A: Math. Theor. 43 (2010), 175303 (13 pages) 



26 



T. Azizov and C. Trunk, "PT symmetric, Hermitian and P-self-adjoint operators related to 
potentials in PT quantum mechanics", J. Math. Phys. 53 (2012), 012109 (18 pages) 

I.V. Barashenkov, S.V. Suchkov, A. A. Sukhorukov, S.V. Dmitriev and Yu.S. Kivshar, 
"Breathers in PT-symmetric optical couplers", Phys. Rev. A 86 (2012) 053809 (12 pages) 

O. Bendix, R. Fleischmann, T. Kottos, and B. Shapiro, "Exponentially fragile PT symmetry 
in lattices with localized eigenmodes", Phys. Rev. Lett. 103 (2009), 030402 (4 pages). 

C. M. Bender, "Making Sense of Non-Hermitian Hamiltonians" , Rep. Prog. Phys. 70 (2007) 
947-1018. 

C. Chong, D.E. Pelinovsky, and G. Schneider, "On the validity of the variational approximation 
in discrete nonlinear Schrodinger equations", Physica D 241 (2012), 115-124. 

S.V. Dmitriev, A.A^ Sukhorukov, and Yu.S. Kivshar, "Binary parity-time-symmetric nonlinear 
lattices with balanced gain and loss". Opt. Lett. 35 (2010), 2976-2978. 

R. Driben and B.A. Malomed, "Stability of solitons in parity-time-symmetric couplers". Opt. 
Lett. 36 (2011), 4323-4325. 

R. Driben and B.A. Malomed, "Stabilization of solitons in PT models with supersymmetry by 
periodic management ", EPL 96 (2011), 51001 (6 pages). 

E. M. Graefe, H. J. Korsch, and A. E. Niederle, Quantum-classical correspondence for a non- 
Hermitian Bose-Hubbard dimer, Phys. Rev. A 82 (2010) 013629 (16 pages). 

M. Hiller, T. Kottos, and A. Ossipov, "Bifurcations in resonance widths of an open Bose- 
Hubbard dimer", Phys. Rev. A 73 (2006) 063625 (4 pages). 

P.G. Kevrekidis, The Discrete Nonlinear Schrodinger Equation: Mathematical Analysis, Nu- 
merical Computations and Physical Perspectives, Springer-Verlag (Heidelberg, 2009). 

O.N. Kirillov, "PT-symmetry, indefinite damping and dissipation-induced instabilities", Phys. 
Lett. A 376 (2012), 1244-1249. 

V.V. Konotop, D.E. Pelinovsky, and D.A. Zezyulin, "Discrete solitons in PT-symmetric lat- 
tices", EPL 100 (2012), 56006 (6 pages). 

K. Li and P.G. Kevrekidis, "PT-symmetric oligomers: analytical solutions, linear stability, and 
nonlinear dynamics", Phys. Rev. E 83 (2011), 066608 (7 pages). 

K. Li, P.G. Kevrekidis, B.A. Malomed, and U. Giinther, "Nonlinear PT-symmetric plaquette", 
J. Phys. A Math. Theor. 45 (2012) 444021 (23 pages) 

A. Mostafazadeh, "Pseudo-Hermitian Representation of Quantum Mechanics," Int. J. Geom. 
Methods Mod. Phys. 7 (2010), 1191-1306. 

Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. N. Christodoulides, Optical Solitons 
in PT Periodic Potentials, Phys. Rev. Lett. 100 (2008) 030402 (4 pages). 



27 



[24] S. Nixon, L. Ge, and J. Yang, "Stability analysis for solitons in PT-symmetric optical lattices", 
Phys. Rev. A 85 (2012), 023822 (10 pages). 

[25] D.E. Pelinovsky, P.G. Kevrekidis, and D. Frantzeskakis, "Stability of discrete solitons in non- 
linear Schrodinger lattices", Physica D 212 (2005), 1-19. 

[26] D.E. Pelinovsky, P.G. Kevrekidis, and D.J. Frantzeskakis, "PT-symmetric lattices with ex- 
tended gain/loss are generically unstable", EPL 101 (2013), 11002 (6 pages). 

[27] D. Pelinovsky and A. Sakovich, "Internal modes of discrete solitons near the anti-continuum 
limit of the dNLS equation", Physica D 240 (2011), 265-281. 

[28] H. Ramezani, T. Kottos, R. El-Ganainy, and D.N. Christodoulides, "Unidirectional nonlinear 
PT-symmetric optical structures", Phys. Rev. A 82 (2010), 043803 (6 pages). 

[29] C.E. Riiter, K.G. Makris, R. El-Ganainy, D.N. Christodouhdes, M.Segev, and D. Kip, "Ob- 
servation of parity-time symmetry in optics", Nature Physics 6 (2010) 192-195. 

[30] S.V. Suchkov, B.A. Malomed, S.V. Dmitriev and Yu.S. Kivshar, "Solitons in a chain of parity- 
time-invariant dimers", Phys. Rev. E 84 (2011), 046609. 

[31] A. A. Sukhorukov, S.V. Dmitriev, S.V. Suchkov, and Yu.S. Kivshar, "Nonlocality in PT- 
symmetric waveguide arrays with gain and loss" Opt. Lett. 37 (2012) 2148-2150. 

[32] A. A. Sukhorukov, Z. Xu, and Yu.S. Kivshar, "Nonlinear suppression of time reversals in PT- 
symmetric optical couplers" Phys. Rev. A 82 (2010), 043818 (5 pages). 

[33] S. Weigert, "Detecting broken PT-symmetry" , J. Phys. A: Math. Gen. 39 (2006), 10239-10246. 

[34] D.A. Zezyulin and V.V. Konotop, "Nonlinear modes in finite-dimensional PT-symmetric sys- 
tems" Phys. Rev. Lett. 108 (2012), 213906 (5 pages). 



28 



